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Abstract. In a transiting planetary system, the presence of a second planet 
will cause the time interval between transits to vary. These transit timing varia- 
tions (TTV) are particularly large near mean-motion resonances and can be used 
to infer the orbital elements of planets with masses that are too small to detect 
by any other means. I present the results of a study of simulated data where I 
show the potential that this planet detection technique has to detect and char- 
acterize secondary planets in transiting systems. These results have important 
ramifications for planetary transit searches since each transiting system presents 
. an opportunity for additional discoveries through a TTV analysis. I present such 

an analysis for 13 transits of the HD 209458 system that were observed with the 
\Q | Hubble Space Telescope. This analysis indicates that a putative companion in a 

. low-order, mean-motion resonance can be no larger than the mass of the Earth 

and constitutes, to date, the most sensitive probe for extrasolar planets that or- 
bit main sequence stars. The presence or absence of small planets in low-order, 
mean-motion resonances has implications for theories of the formation and evo- 
lution of planetary systems. Since TTV is most sensitive in these regimes, it 
\ should prove a valuable tool not only for the detection of additional planets in 

transiting systems, but also as a way to determine the dominant mechanisms of 
planet formation and the evolution of planetary systems. 

> 

1. Introduction 

The core accretion model of planet formation often predicts the presence of 
small planet masses that become trapped in mean- motion resonance with a gas 
giant planet as the giant planet migra t es inward (Tnommes|l2005l ; Zhou et al 
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Or 



20051 : 1 Papaloizou fc Szuszkiewiczl 120051 ; iRavmond et al.ll2006l ) . The small sizes 



and masses of these putative trapped planets render them difficult to detect 
by conventional methods like radial velocity (KV) mea sure ments and planetary 
transi ts. However, recent work by lAgol et ahl ( 20051 ) and Holman &: Murray! 



(2005) indicate that, in a known transiting system, the times of the planetary 
transits are a particularly sensitive probe for additional, resonant bodies. 

If a second planet exists in a transiting system, then dynamical interactions 
within the system cause the time interval between transits to vary. For systems 
in mean-motion resonance, these transit timing variations (TTV) can be quite 
large; often of order (m^/m^P^ in where the m p and tot are the masses of 
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the perturbing and transiting planets respectively and is the period of the 
transiting planet (this amounts to nearly 15 minutes for an Earth-mass perturber 
and a Jupiter-mass transiting planet that is in a 3-day orbit )Q. From this TTV 
signal, one may infer the orbital elements and mass of the perturbing planet. 

The strong sensitivity of TTV to resonant systems should allow for the dis- 
covery of sub Earth-mass planets in transiting systems, possibly using modest 
(~ 1 m) ground-based telescopes. Indeed, the first TTV analysis was conducted 
for the TrES-1 system where ground-based observations of 1 1 transits were 
able to probe for planets with masses smaller than the Earth (jSteffen fc Ago! 
20051 1. Here we present portions of a similar study of Hubbl e Space Telescope 



(HST) observations of 13 transits of the HD 209458 system (jBrown et alj|2001 



Schultz et afll2004 iKnutson et ai1l2006l ) . We compare the sensitivity to mass of 
ground-based observations with space-based observations. We also present the 
results of an analysis of 100 simulated systems where we attempt to infer the 
orbital elements of each system from a set of transit times. 



2. Analysis of HD 209458 System 

We analyzed 13 transit times that were obtained from observations of HD 209458 
using various instruments on HST. The details of the photomet ric data reduc- 
tion, t he transit times, and the balance of our results appear in lAgol fc Steffenl 



(120061 ) . While those data give no evidence for additional planets, we use them 



to constrain the presence of secondary planets as a function of various orbital 
elements. Of particular interest are the constraints on planets in the 2:1 mean- 
motion resonances (see Figure [1]). We find that a secondary planet in the exterior 
resonance cannot be much more massive than the Earth. For the interior reso- 
nance, a second planet must be significantly less massive than the Earth in order 
for the systems to be consistent with the data at the 3-<r level. 

In addition, we simultane ously analyzed t hese t ransit data and 68 of the 
RV measurements reported by lLaughlin et aL ( 20051 ) in effort to constrain the 



allowed mass of a second planet. Here, we restrict our analysis to initially 
circular orbits. Figure [5] shows the results of this an alysi s and compares it wit h 
the theoretical predictions given in lAgol et al.1 420051 ) and lSteffen fc Aeoll (120051 ). 



We find that for many orbits, particularly those where the period ratio is larger 
than 2:1, the RV technique is more sensitive than TTV. It also illustrates the 
tremendous advantage that TTV has in probing resonant systems. 

Our TT V analysis cannot exclude the putative secondary planet that was 
proposed by iBodenheimer et all ( 20031 ) as an explanation for the apparently 



large radius of HD 209458b. The proposed planet is relatively distant (having a 
period of ~ 80 days) and is not near a mean-motion resonance. Consequently, 
it would induce a TTV signal with an amplitude less than one second. RV 
measurements are better suited to test this hypothesis which we did not attempt. 



1 We note that lAgol et~afl (|2005T ) and lHolman fc Murray! lj2005l ) define the TTV signal differently. 
The former being the residuals from a linear fit to the period and epoch, the latter being the 
difference between successive periods. The former gives a larger signal for resonant systems 
while for non-resonant systems the signals are comparable. When analyzing data, one will 
likely fit the transit times directly and will thus render this difference immaterial. 
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Figure 1.: Maximum allowed mass of a secondary planet in the HD 209458 system 
near the 2:1 interior (upper panel) and exterior (lower panel) resonances for any 
value of eccentricity. The dashed lines correspond to the mass of the Earth. 




Figure 2.: Maximum allowed mass for a second planet on an initially circular 
orbit in the HD 209458 system as a function of its period ratio with the known 
plan et. The gray cu r ve is the result from a TTV analysis of the transit times 
from lAgol fc Steffenl (120061 ) the dotted curve is the theoretical prediction from 
the appendix of Agol et al. ( 20051 ). the black curve is the result of simultaneousl y 
analyzing the transit times and the RV measurements of lLaughlin et al.l ((2005), 
the t hick dashed line is th e prediction for the sensitivity of RV me asurements 
fromlSteffen & Agoll (|2005l ). the large dots are given by equation 33 of I Agol et al.l 
(|2005F and the horizontal dashed line corresponds to the mass of the Earth. 



3. Ground-based vs. Space-based 



The tr ansit times for the TrES-1 system that were published bv lCharbonneau et al 
(2005) have an average timing error of about 100 seconds. The most precise of 
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the times, derived from observations with the 1.2m Fred Whipple telescope, had 
an uncertainty of 26 seconds. In contrast, the average uncertainty for the HST 
observations of HD 209458 is 25 seconds with the best being 10 seconds. This 
serves to illustrate the opportunities available to ground-based facilities. If the 
timing precision of 26 seconds is reproducible, then a set of such observations 
could readily probe fo r resonant planets with masses smaller than the Earth. 

Recent studies bv lHolman fc Winnl ( 20061 ) indicate that ~ 10 second timing 
precision should be possible with ground-based telescopes. This fact, and the rel- 
ative ease of obtaining time on such telescopes compared with their space-based 
counterparts, shows that ground-based, follow-up observations of transiting sys- 
tems identified by Kepler and CoRoT should prove invaluable. Those observa- 
tions could address important issues in planet formation and evolution and may 
identify the smallest known extrasolar planets for the foreseeable future. 



4. Identifying Planets with TTV 

Once a TTV signal is identified there remains the question of how difficult it 
is to infer the orbital elements of the perturbing planet. We outline here our 
preliminary investigation of this issue. For this study we randomly generated 
100 systems that are near the 2:1 mean- motion resonance where the perturbing 
planet is exterior to the transiting planet. The distributions from which we drew 
our orbital elements are uniform for the longitude of pericenter, time of pericen- 
ter passage, and period ratio (within 5% of the resonance) and logarithmically 
uniform in eccentricity and mass. For the unknown planet, the eccentricities are 
between 10 -3 and 1 and the masses are between 10 _6 Mo and 10 _3 Mo where 
Mo is the mass of the central star. The mass of the known planet was fixed at 
10~ 3 Mo and its eccentricities are between 10~ 3 ' 5 and 10~ ' 3 . We rejected any 
system where the hill spheres of the two planets could touch. We chose resonant 
systems because that is the regime where the TTV technique is unique in its 
sensitivity and where it will likely prove most valuable. We allow the perturbing 
planet masses to be so large (large enough to be seen with RV measurements) 
because it may be possible to conduct a TTV analysis of Kepler and CoRoT 
data before RV measurements are available. 

For each system we generate a set of 120 consecutive transit times. To each 
time we add Gaussian, white timing noise with a dispersion equal to 5 seconds 
per day of the period of the transiting planet (e.g. for a 3-day period transiting 
planet the timing error is 15 seconds). From these transit times we analyze 5 
subsets of data: (1) all 120 transit times, (2) 60 randomly selected times, (3) 
30 random times, (4) 15 random times, and (5) 15 correlated times where, if a 
particular transit is selected then there is a 67% chance that a second transit is 
selected within 4 orbits and, if a second transit is selected, then there is a 50% 
chance that a third transit is selected within 4 orbits of the second transit. This 
correlated selection should roughly approximate the effects of realistic programs 
which are constrained by observation conditions, seasonal variations, etc. though 
these criteria were not selected with any serious rigor. 

For our analysis, we assumed that the orbital elements of the transiting 
planet are known and that the orbits of the two planets are coplanar. We sys- 
tematically step through the period ratio of the planets between the 3:2 and the 
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Correct Identifications in 2:1 Exterior Resonance 
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Figure 3.: The fraction of correctly identified planets in the exterior 2:1 mean- 
motion resonance with a known planet for analyzes of 120, 60, 30, 15, and 15 
correlated transit observations. The bar on the left of each set is for planets 
with masses less than the mass of the Earth. The center bar of each set is for 
planets with masses between one and ten Earth masses. The bar on the right is 
for planets larger than 10 times the mass of the Earth. 



7:1 resonances and then marginalize over the eccentricity, longitude of percen- 
ter, time of pericenter passage, and mass of the perturber at each point. This 
marginalization is done by choosing 500 random values for the parameters at 
each value of the period ratio, generating hypothetical transit times for each sys- 
tem, calculating the x 2 between the actual transit times and the hypothetical 
transit times, and requiring the set of parameters with the lowest \ 2 to represent 
the best-fit system for that value of the period ratio. 

For each system and each set of transit timing data we consider the iden- 
tification of the perturbing planet to be successful if one of the three lowest 
values of the x 2 — that is, the x 2 as a function of period ratio — is within 5% of 
the actual period ratio. Often the remaining orbital elements are also correctly 
identified, but for our purposes we neglect this information since, once the period 
ratio is identified a more detailed analysis could be done to find the remaining 
parameters and their uncertainties. Our results are shown in Figure [3] where we 
see that a significant fraction of the systems are correctly identified. 

These preliminary results likely represent a lower bound for the fraction 
of systems that can be identified from such data. A more sophisticated search 
algorithm will certainly do better than the one that we use, which is chosen for 
simplicity and speed and whose scope is only to provide some basic insight into 
the challenges associated with inferring the parameters. Some common charac- 
teristics of the systems which we fail to identify are that the perturber masses 
are often small, the eccentricities are often large, and the period ratios that have 
the lowest x 2 are correspond to resonances other than the 2:1 resonance. Such 
resonances often have similar TTV signals over the time span of the data. 

We mention one issue exposed by this experiment. For a TTV analysis of 
transit data it is unlikely that a basic Markov Chain Monte Carlo (MCMC) 
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algorithm will be appropriate. This is because, for a system in resonance (where 
TTV is most useful) the x 2 often has many deep, narrow local minima that 
correspond to various resonances. Such minima are difficult for a Markov Chain 
to handle because, if the step size of the chain is too small it will become trapped 
in these local minima. Yet, if it is large enough to escape the local minima then 
it is often too large to properly resolve them. Thus, a multi-stage investigation 
will probably prove more fruitful where, for example, the local minima are first 
identified using a search similar to the one outlined above and then each minima 
is studied with MCMC such that the chain remains inside the selected minima. 

5. Discussion 

The results presented in this article indicate that TTV should be a valuable 
tool for both discovering additional planets and for constraining the relative 
importance of different planet formation and evolution mechanisms. A set of 
a few tens of high quality, ground-based transit observations can readily probe 
for Earth-mass or smaller planets in several mean- motion resonances. A similar 
number of transit observations with space-based telescopes should be able to 
identify mars-mass companions. Such discoveries could place constraints on the 
initial mass function of the planetesimals and would provide important guidance 
to theorists who model planet formation. We showed that relatively simple 
search algorithms can correctly identify a large fraction of possible resonant 
systems. Thus, if small resonant bodies do exist in planetary systems we are 
confident that they can be found using the TTV technique. 
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